APS/123-QED 



Soft gamma-ray background and light Dark Matter annihilation 

Yann Rasera, Romain TeyssierQ and Patrick Sizim 
CEA Saclay, Bat. 709, 91191 Gif-sur-Yvette, France 

Michel Casse 

CEA Saclay, Bat. 709, 91191 Gif-sur-Yvette, France and 
Institut d'Astrophysique de Paris, 98bis Bd Arago, 75014 Paris, France 

Pierre Fayet 

Laboratoire de Physique Theorique de I'ENS, UMR 8549 CNRS, 
24 rue Lhomond, 75231 Paris Cedex 05, France 

Bertrand Cordier and Jacques Paul 
CEA Saclay, Bat. 709, 91191 Gif-sur-Yvette, France 
(Dated: March 9, 2006) 

The bulk of the extragalactic background between 10 keV and 10 GeV is likely to be explained 
by the emission of Seyfert galaxies, type la supernovae, and blazars. However, as revealed by 
the INTEGRAL satellite, the bulge of our galaxy is an intense source of a 511 keV gamma-ray 
line, indicating the production of a large number of positrons that annihilate. The origin of the 
latter is debated, and they could be produced, in particular, by the {S- or P-wave) annihilation 
of light Dark Matter particles into e'*'e~. In any case, the cumulated effect of similar sources at 
all redshifts could lead to a new background of hard X-ray and soft gamma-ray photons. On the 
basis of the hierarchical model of galaxy formation, we compute analytically the SNIa contribution 
to the background, and add it to Seyfert and blazars emission models. Confronting these expected 
contributions to observation, we find that any extra contribution to this unresolved background 
around 511 keV should be lower than about 4 keV cm^^ s^^ sr^^. We also estimate analytically the 
extragalactic background due to Dark Matter annihilation, increasing the accuracy of the earlier 
computations. Indeed, we take into account the large positron escape fraction from low mass dark 
matter halos, unable to confine a dense and magnetized interstellar medium. Our new background 
estimate turns out to be one order of magnitude lower so that the hypothesis of a light Dark 
Matter candidate remains compatible with the observed extragalactic background for a wider range 
of particle masses and cross-sections. 

PACS numbers: 95.35.-|-d, 95.85.Nv, 95.85.Pw, 97.60.Bw 



I. INTRODUCTION 

The cosmic gamma-ray background (CGB) between 10 
keV and 10 GeV has been measured by several gamma- 
ray satellites (HEAO, SMM, COMPTEL and EGRET) 
Below 100 keV, it is believed that the main con- 
tribution comes from Seyfert galaxies |1|| j48||. Above 10 
MeV, a simple model for blazars reproduces both the 
amplitude and the slope of the data In the interme- 
diate energy range, however, another type of sources is 
needed, since blazar spectra show a clear break near 10 
MeV and the cosmological gamma-ray background from 
Seyfert galaxies falls off above about 100 keV (see Fig.^). 
As discussed by several authors "T] , type la super- 
novae could make a significant contribution in this energy 
range, which we shall evaluate in Section 2. 

Furthermore, the recent observation, by the INTE- 
GRAL satellite, of a 511 keV diffuse emission line from 
the galactic bulge Q shows that electron-positron anni- 
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hilation are taking place there with a very large rate 
~ 1.5 10'*'^ s~^. Such a high rate is difficult, if not impos- 
sible, to impute to astrophysical objects, and the source 
of these positrons in the bulge is subject to intense de- 
bate. This emission from the bulge of the Milky Way 
could be the signature of light Dark Matter particles an- 
nihilating into e+e~ 0,0, the positrons eventually an- 
nihilating with electrons encountered in the interstellar 
medium ,11.] . In any case, if one extrapolates this diffuse 
emission to all other galaxies in the Universe, the inte- 
grated flux could make a significant additional contribu- 
tion to the hard A"-ray and soft gamma-ray background. 

Since the potential implications of such an hypothesis 
for particle physics and cosmology are very important, we 
want to test here its validity or at least its consistency, 
as far as the cosmological gamma-ray background is con- 
cerned. Using a recent model of galaxy formation '12], 
we compute self-consistently the gamma-ray background 
coming from both type-la supernovae (adding Seyfert- 
galaxies and blazars), and annihilating positrons from 
light Dark Matter halos. We follow in this respect the 
earlier work of Ahn & Komatsu (2| , who were the first 
to compute this background, assuming, like them, how- 
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FIG. 1: Diffuse background spectrum as a function of photon 
energy inspired by Fig. 4 of Strigari et al. 6] (we remove the 
error bars for sake of visibility). The crosses (HEAO), stars 
(COMPTEL) and diamonds (EGRET) correspond to the ob- 
servations Q, Q . At low energy, Seyfert galaxies (dashed line) 
are the main contributors I'll. At intermediate energy, Type 
la supernovae (continuous line), as calculated in this article, 
dominate. At high energy, blazars (dot-dashed line) explain 
the observed cosmological gamma-ray background JJ. Alto- 
gether the sum of the three contributions (dotted line) is a 
factor of 2 below the observations near 511 keV (indicated 
by an arrow). An additional contribution from light Dark 
Matter particles, of up to about 4 keV cm~^ s~^ sr~^, is not 
excluded. 



ever, that positrons annihilate "on the spot" in all Dark 
Matter halos in the Universe. These authors explored 
various scenarios for the internal structure of Dark Mat- 
ter halos and analyzed how these various models affect 
the amplitude of the gamma-ray background. 

In the present paper, we would like to go one step fur- 
ther, and explicitly take into account the role of baryons 
in the process of positron confinement and annihilation. 
This process is indeed possible only if the parent halo 
contains enough baryons (and therefore also electrons) 
to host a dense, magnetized, interstellar medium. This 
is a necessary condition for the "on the spot" approxima- 
tion to be valid. Since low mass halos are unable to host 
enough baryons in a cold and magnetized disc [T^. the 
escaping positron mean free path increases dramatically. 
To compute the positron escape fraction and their propa- 
gation in the expanding background, a complex diffusion 
study would be necessary. As a first order approxima- 
tion we shall assume that the positron escape fraction 
goes from zero to one below the critical mass for a Dark 
Matter halo to host a galaxy, as computed in and 
that these escaping positrons never annihilate. 

The outline of this article is as follows. In Section 2, 



we estimate the contribution of SNIa to the gamma-ray 
background, adding it to that of Seyfert galaxies and 
blazars (at lower and higher energies, respectively). The 
difference with the observed spectrum provides an up- 
per limit on any additional contribution such as the one 
due to annihilating positrons, that could come from light 
Dark Matter particle annihilation. In Section 3, we cal- 
culate the diffuse cosmological background induced by all 
Dark Matter halos in the Universe, taking into account 
that positrons cannot annihilate in small mass halos, and 
compare it to the previous calculation performed by Ahn 
& Komatsu p^. In Section 4, we present the gamma-ray 
background constraints on the annihilation cross-sections 
and the masses of light Dark Matter candidates. We 
summarize our main conclusions in Section 5. In an Ap- 
pendix, we estimate for calibration purpose, and confront 
with SPI/INTEGRAL observations, the 511 keV emis- 
sion from the galactic bulge. We consider various annihi- 
lation cross-sections (depending on whether they are S- 
or P- wave dominated) and the corresponding halo Dark 
Matter density profiles. 

II. DIFFUSE GAMMA-RAY BACKGROUND 
FROM SNIA 

The type la supernovae contribution to the gamma- 
ray background depends primarily on the star formation 
history in the Universe, which, in this paper, is derived 
from a new self-consistent model of galaxy formation 12| . 
This analytical model predicts the cosmological evolu- 
tion of the four main baryon phases in the Universe: dif- 
fuse intergalactic gas, hot gas, cold gas in galaxies and 
stars. These theoretical predictions were validated with 
high-resolution cosmological simulations using the RAM- 
SES [l| and GADGET codes They also reproduce 
the observed amount of cold gas in the Damped Lyman- 
Alpha systems_T&jj|_^J_and the observed Cosmic Star For- 
mation Rate lilM M m m m El m (see Fig. Ell, 
which is of prime interest here. The model particularly 
emphasizes the important cosmological role of the mini- 
mal mass for a halo to host galaxies, M^niniz). 

We shall use it to, first, compute the SNIa gamma- 
ray background, and, also, to evaluate the new back- 
ground that could be attributed to the annihilation of 
positrons, possibly generated in annihilating Dark Mat- 
ter halos. This internal consistency allows us to perform 
a fair comparison between the two types of gamma-ray 
sources. 

In the general case, the background intensity is given 
by 

h = j,(//(l + z), z) dt , (1) 

with j,j{i',z) the comoving emissivity at redshift z, or 
time t. 

The gamma-ray lines from SNIa result from the explo- 
sive synthesis of radioactive ^^Ni nuclei, decaying succes- 
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FIG. 2: Comoving supernova rates as a function of redshift 
from our analytical model. The dashed line shows the SNII 
rates and corresponds to the star formation rate multiplied by 
the fraction of SNII per unit of stellar mass formed {esNii — 
0.007 Mg^). The dot-dashed line represents the SNIa rate 
used for computing the diffuse gamma-ray background. Both 
rates are compatible with the observations from Dahlen et al. 

f (symbols) and with the upper limit from Strigari et al. 
(dashed region). 

sively into '^^Co and ^^Fe. The comoving emissivity from 
SNIa may be expressed as 



~ hv p^{t-tsN) ^SN 
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with (t — Isn) the comoving star formation rate at 
time t — tsN from our analytical model (see Fig. [21), tsN 
denoting the average delay between star formation and 
SNIa explosion, csn is the number of SNIa per unit of 
stellar mass formed, M^j the mass of Nickel produced and 
'S'sAr(j^) the average spectrum per Nickel nucleus. We 
take tsN ^ 2.5 Gyr, esN = 1-4 x lO'^ Mg^ (so that 
the resulting SNIa rates be within the 2a error bars of 
observed rates H^), Mp , ~ 0.5 Mq, and the spectrum 
S{i') as computed in [2g. 

The resulting extragalactic background spectrum from 
SNIa is presented in Fig. ^ It shows a bump in the 
range 300 keV up to 3 MeV, at a level which turns out 
to be close to the predicted contributions from Seyfert 
galaxies and blazars. This contribution, in agreement 
with the SNIa contribution from Strigari et al. Q, is 
slightly higher than their preferred model ji^, because 
our star formation history is sli ghtl y more efficient, as 
suggested by recent observations [ij . 

Altogether the resulting evaluation from known as- 
trophysical sources reproduces reasonably well the ob- 
served extragalactic background below 100 keV and 



above 3 MeV. However, in the range from 100 keV to 
3 MeV the three contributions fall short of explaining the 
bulk of the Cosmic Gamma-Ray background (as empha- 
sized by Strigari et al. Particularly, in the 100 keV- 
511 kcV range of interest in this article, the sum of the 
three contributions appears to be lower with a difference 
of the order of 4 keV cm~^ s~^ sr~^. This sets an upper 
limit on a possible Dark Matter annihilation signal. 



III. DIFFUSE BACKGROUND FROM 
COSMOLOGICAL HALOS 

A. Diffuse background 

The diffuse background is simply the sum of the red- 
shifted emissions from positron annihilation in all cos- 
mological halos, in principle at all redshifts (Eq. The 
comoving emissivity can be computed by summing up 
the individual halo emissivities 
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where Mmin is the minimal mass for emitting halos. The 
luminosity per halo integrated up to the radius of the 
halo i?2oo is 



Pu (r) 47r dr 



(4) 



and N{M, z) is the Press-Schechter .27| mass function 
for cosmological halos. Considering positron annihilation 
"on the spot" , the volume emissivity is given by 



(5) 



'■X 



with px{r) the Dark Matter mass density profile, 
(f Wrci ('')) is the annihilation cross-section, the factor i 
being present only in the case of self-conjugate Dark Mat- 
ter particles [50l |. 

Positronium annihilation introduces a specific emission 
spectrum S'pos(i'), with 25% of the energy injected in the 
511 keV line, and the remaining 75% spread over a 37 
continuum. 



B. Dark matter density profile 

The mass distribution in each Dark Matter halo is in 
fact quite uncertain. Ahn and Komatsu T^l have ex- 
plored a wide range of halo density profile parameters. 
In this paper, we restrict ourselves to Dark Matter dis- 
tribution parameters at face value, as suggested by TV- 
body simulations, based on the following general fitting 
formula 



Px{r) (X 



x'y{l + x") 



(6) 
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where x = rjrs with the scahng radius corresponding 
to the concentration parameter c = R2m/fs (typically 
between 4 and 40 depending on halo mass and redshift). 
7, a and [3 control the slope respectively for small (r < 
Ts), intermediate (r ~ r^) and large radii (r > r^). 

The concentration parameter defines whether halos are 
rather peaked (c ~ 40) or shallow (c ~ 4). Here again 
the mean value as a function of redshift and halo mass is 
given by a fit on cosmological simulations [28|. 

c = maa;(4,4ii^) (7) 
1 + z 

with Zc the collapse redshift given by M*(zc) = lO^^M, 
A/* being the non- linear mass at a redshift z. This for- 
mula is valid only for halos greater than ~ 10^ M0. The 
behaviour of the concentration parameter for smaller halo 
masses, unresolved by numerical simulation, is totally un- 
known. 

The slope at the center of Dark Matter halos is thought 
to be between a = 1 and a = 1.5. We therefore consider 
two extreme dark matter profiles as given by Navarro, 
Frenk and White "291 (a 1, /3 = 3 and 7 = 1) and 
Moore ISflj {a = 1.5, /? = 3 and 7 = 1.5). Note how- 
ever that both density profiles saturate at very low ra- 
dius Rmin when n{Rrmn) (crurei (r)) = 1/iff due to self- 
annihilation {n{Rmin) = Px{Rmin) /fTix is the dark mat- 
ter numeric density and tn is the age of the universe). 



C. Annihilation cross-section 

The relic abundance of Dark Matter particles depends 
on their decoupling temperature, which is a fraction of 
their mass (i.e. Tp = with xp ~ 16 to 20 depending 
on mx) and is, roughly, inversely proportional to their 
annihilation cross-section {avrc\/c)p at freeze-out. The 
values required for a correct abundance, corresponding 
to ildm — 23 %, are then, for such light particles, of 
the order of a few (up to « 10) picobarns (correspond- 
ing to (crwrci)i? ~ 10~^^ cm^ s~^) 0, 0, depending 
on whether they are self-conjugate or not, and on the 
possible velocity-dependence of their annihilation cross- 
section at freeze-out. 

Such values are in any case rather large compared to 
ordinary weak-interaction cross-sections, especially when 
dealing with light particles. This necessitates an unusual, 
more powerful, annihilation mechanism, that could result 
from the exchanges of a new light neutral gauge boson 
U , or, in the case of spin-0 Dark Matter particles^from 
the exchanges of new heavy (e.g. mirror) fcrmions |9lll(l 

mm 

A rather large annihilation cross-section could lead to 
an excessive continuum of gamma-ray photons at var- 
ious energies (depending on mx)- Cross-sections which 
behave, at least to a large extent, proportionally to v"^ [P- 
wave annihilation), may therefore be preferred 32], es- 
pecially at lower mx ■ The residual annihilation of Dark 



Matter particles in bulges of spiral galaxies or in ellipti- 
cals would then include a suppression factor that could 
be, in the pure P-wave case, as strong as Whaio/^F ~ 
10^^. Furthermore, and independently of the above ar- 
gument, lighter Dark Matter masses mx tend to be pre- 
ferred, to avoid excessive gamma-ray production as com- 
pared to production, in our galaxy |lfl l33||. 

We shall therefore consider annihilation cross-sections 
parametrized as crwroi — a+hv^ ,vi\th. {a v^c\) p ~ 10~^^ 
cm"^ s~^ at freeze-out (for a self-conjugate particle - or 
twice this value, for a non self-conjugate one). And ex- 
plore in particular, for low- velocity halo particles, the two 
extreme situations (crvroi) — a (S'-wave) and ~ hv^ (P- 
wave annihilation). The resulting emission profiles de- 
duced from a given Dark Matter profile px are computed 
from (awrei) Pxl'^x (^f- Eq. El for a self-conjugate par- 
ticle [53 )• See also the Appendix for further comments. 

Note that for the pure P-wave cross-section, the emis- 
sivity now depends on the Dark Matter 3D velocity dis- 
persion We therefore compute a\j-, as a function of 
the radius by solving the Jeans equation for a NFW or 
Moore potential. The resulting emission profiles turns 
out to be less peaked than for a pure 5-wave cross- 
section. 



D. Role of the baryons 
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FIG. 3; Redshift evolution of the minimal halo mass Mmin be- 
low which cold disk gas cannot form (see for details). We 
use a reionization redshift z,. ~ 20 as suggested by WMAP. 

This rather standard approach has been ap plie d to 
compute the soft gamma-ray background in |l3|. in- 
tegrating individual halo emissivity over the Press & 
Schechter distribution (given by Eq. |3Jl, using as lower 
bound of the integration interval the maximum between 
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the Dark Matter free-streaming and the Dark Matter 
Jeans masses. This leads to a minimal mass (M„iin) equal 
to a fraction of solar mass. As a consequence the comov- 
ing gamma-ray emissivity is dominated by the cumulated 
emission of numerous small mass halos (if one uses the 
concentration parameters given by Eq. 13 52] ) . 

However, as we have already discussed, small mass 
halos are unable to retain gas and annihilate positrons, 
and therefore cannot contribute to the gamma-ray back- 
ground in which we are interested. We take into ac- 
count the crucial role of the baryons (and associated 
electrons, and magnetic fields) in confining and anni- 
hilating the Dark Matter positrons. This trapping can 
be achieved only if both the density of the interstel- 
lar medium and the galactic magnetic field are suffi- 
cient. Indeed cosmological simulations of galaxy forma- 
tion in the hierarchical framework of structure formation 
show that baryons cannot collapse and form high-density 
centrifugally-supported gas discs in halos having a mass 
lower than a minimal value Afmin ~ 10^ — 10^^ h^^ Mq. 

This critical mass threshold is a key ingredient of the 
current galaxy formation theory. Gnedin showed that the 
fraction of baryons decreases str ong ly in halos smaller 
than the so-called filtering mass '35'| , as a consequence 
of the non-zero temperature of the intergalactic medium 
which prevents gas from collapsing into too small Dark 
Matter halos. Hoeft et al. [S^ also showed that the halo 
mass must be greater than the minimal cooling mass; if 
not, the fraction of baryons is high but galaxies cannot 
form because cooling is inefficient. The resulting minimal 
halo mass Afmin for galaxy formation is then the maxi- 
mum between the minimal cooling mass and the filtering 
mass, as computed in p^. The evolution of this minimal 
mass with redshift is shown in Fig.|3| It is of course much 
larger than the one used in TS] . 

As computing accurately the escape fraction of 
positrons as a function of halo mass is beyond the scope 
of this paper, we shall consider here, for simplicity, that 
below Afmin, essentially all positrons escape the halos 
(and we neglect their contribution to the background), 
while above Afmin, confinement is supposed to be effi- 
cient and all positrons are taken to contribute. Further- 
more, the contribution from diffuse baryons in the Uni- 
verse (baryons which are not in collapsed halos) is neg- 
ligible in this range of wavelength because the annihila- 
tion time scale is larger than the age of the Universe ■ 
Note that, while the "on the spot" approximation, used 
by Ahn & Komatsu for the whole mass range, leads, 
in our opinion, to an overestimation of the diffuse back- 
ground, our approach, though more accurate, should lead 
to an underestimation of the background level. 

Fig.^lillustrates the distribution of comoving positron 
emissivity as a function of halo masses for different red- 
shifts. The comoving emissivity is dominated by small 
mass halos. However, only a small fraction of the halos 
contains enough gas and magnetic field for the produc- 
tion of gamma-rays: altogether, only a fraction of about 
10 % of the emitted positrons are converted into photons. 
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FIG. 4: Distribution of comoving positron emissivity as a 
function of halo masses for three different redshifts: z — 
(continuous line), z = 6 (dashed line) and z = 12 (dot-dashed 
line). For each redshift, the vertical line indicates the mini- 
mal mass for a halo to host a galaxy and therefore to be able 
to annihilate positrons into gamma-rays. The bold lines de- 
lineate the regions, to the right of the vertical lines, where 
positrons are converted into photons. It follows that only 
about 10 % of the released positrons are actually converted 
into gamma-rays. 



As a consequence, our computation of the gamma-ray 
background leads to a result about 10 times smaller than 
the one evaluated by Ahn and Komatsu . 

More precisely, we have computed the soft gamma-ray 
background for direct annihilation with a S-wave cross- 
section, a NFW Dark Matter density profile and a parti- 
cle masse mx = 20 MeV in the two different cases : dark 
matter-based minimal mass as in [T^ and baryon-based 
minimal mass as in the present paper. Fig. |S1 illustrates 
the decrease by a factor of 10 using our new approach. 
Thus, if we would like to reach the same level of back- 
ground, we would have to divide the particle mass mx 
by a factor of ~ \/T0. 

Note that the spectral shape is also modified. The 
spectrum declines at low energy because it corresponds 
to high redshift where halos greater than the minimal 
mass for galaxy formation becomes rare. Note that in 
Ahn and Komatsu 13], identical results were recovered 
for the same Dark Matter particle mass, using rather ex- 
treme values for the halo concentration parameter c. Re- 
call that here we consider halo concentration parameters 
only at face value, as predicted by iV-body simulations, 
but we take into account baryon physics as part of the 
annihilation mechanism. 
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FIG. 5: Cosmic gamma ray spectrum produced by NFW pro- 
files witii mx = 20 MeV, a S-wave cross-section and, consid- 
ering that no positronium is formed at all (for the purpose of 
comparison with earlier results) . The upper dot-dashed curve 
was computed using the dark matter-based minimal mass (as 
in Ahn and Komatsu p^l. On the contrary, for the lower 
continuous curve we took into account baryon physics, intro- 
ducing a minimal mass for a halo to host a galaxy. The crosses 
and stars are the observational data. 

E. Calibration on Milky Way 

As we have seen in the previous sections, the diffuse 
gamma-ray background depends on three main quanti- 
ties. The first is the annihilation cross-section: we are 
going to explore two extreme cases: S'-wave and P-wave. 
The second ingredient is the dark matter mass density 
profile: we are going to test peaked distributions (Moore, 
c = 15) and shallow ones (NFW, c = 5). The last un- 
known quantity is the dark matter particle mass mx- 
In this Section, we are going to fix the latter (for a given 
set of annihilation cross-section and density profile) using 
the constraint set by the detected galactic signal. 

The fine emission at 511 keV detected by INTEGRAL 
from the galactic center region is at a level of 10^'^ ph 
cm~^ s~^ Q. Several types of astrophysical sources have 
been considered as potential candidates to explain this 
emission. However, SNIa fall short sustaining the high 
positron injection rate '37, "s^. Hypernoyae [s^ and the 
related gamma-ray bursts J£, JQ, ,4l| are in a better 
position, but since the number of massive stars is about 
ten times larger in the disk than in the bulge, hot spots of 
511 keV emission should show up in the disk plane [s^l, 
which is not the case. Low mass X-ray binaries have 
also been suggested js^, but no 511 keV emission has 
been observed from these objects. We therefore consider 
the hypothesis that light Dark Matter particles annihi- 
late into electron-positron pairs, mainly in the galactic 



mx (MeV) 


S-Wave 


P-wave 




c=15 


1500 


2.4 


Moore 


c=10 


900 


1.2 




c= 5 


440 


0.44 




c=15 


190 


0.42 


NFW 


c=10 


110 


0.20 




c= 5 


45 


0.060 



TABLE I: This Table summarizes, for different cases, which 
dark matter particle mass (in MeV) is (or would be) required 
to reproduce the level of the INTEGRAL signal from the 
galactic bulge. We explore different cross-sections (S-wave 
or P-wave), different inner slopes (Moore or NFW) and dif- 
ferent concentration parameters (c = 5 — 10 — 15). (A result 
smaller than ^ MeV indicates that the signal cannot be repro- 
duced with the annihilation cross-section and density profile 
considered.) 



center region where the Dark Matter density is at a max- 
imum 0, lifl yJI • The resulting low-energy positrons 
are confined by the magnetic field of the bulge, where 
they are progressively slowed down by ionisation losses. 
A large fraction (0.93) forms positronium with ambient 
electrons and annihilate into two (25% of probability) or 
three 7 photons '42*1 . Positronium formation plays there- 
fore an important role because it decreases by a factor of 
about 3 the intensity of the 511 keV line. 

The flux of gamma rays (from a direction making an 
angle 9 with the direction of the galactic center) is given 
by the integral of the emissivity along the line of sight 

FAO) = ^ I Mr) dl , (8) 

where (r) is the volume emissivity. We assume that es- 
sentially no annihilation can take place outside the stellar 
bulge, due to a lack of gas (corresponding roughly to an 
angle 6'buigc = 16°). 

The resulting profiles are convolved with the INTE- 
GRAL/SPI Point Spread Function (PSF). This method 
allows a fair comparison between observations and mod- 
els, and depends only weakly of the poorly known size 
of the gaseous bulge. For each couple of cross-section - 
dark matter profile, we have computed which mass uix 
fits best the observed level of 511 keV emission. The 
results are summarized in Tabled 

As expected, the most peaked profiles (such as Moore 
profile) and the most concentrated ones (c = 15) require 
the largest values of mx , since for a given mass density 
the number density, and hence the annihilation rate de- 
creases when the Dark Matter particle mass increases. 
It is worth noting that there is a factor of 10^ — 10'^ 
for the mass (corresponding to 10^ — 10^ in flux) in fa- 
vor of the S'-wave cross-section compared to the P-wave 
one. Both the NFW S'-wave case and the Moore P-wave 
reproduce the total flux of the bulge 511 keV emission 
with reasonable Dark Matter particle mass of the order 
of mx ~ 100 MeV and mx — 1 MeV, respectively. On 
the opposite, the NFW P-wave case would require masses 
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{mx < 0.42 MeV) so small that they are unable to pro- 
duce 511 keV photons. And the Moore 5- wave case re- 
quired such high masses {mx > 440 MeV) that they 
would lead to an excessive bremsstrahlung emission of 
soft gamma-ray photons js^ . 

Among the specific cases considered two models (NFW 
with S-wetve and Moore with P-wave) are therefore fa- 
vored by the galactic signal. If one considers only the 
total emission from the bulge, neither the S'-wave nor 
the P-wave cross-section can be excluded. If one consid- 
ers the emission profile however, a different conclusion 
could be drawn: according to a recent article [2| the 
shape of the emission profile could be used to exclude the 
P-wave scenario. We address this interesting question in 
the Appendix, in which we conclude that, to our opinion, 
both scenarios cannot be discriminated yet. 



F. Results 



As shown in Fig. the background predictions in the 
two models are at the same level with a slightly differ- 
ent spectral signature. The main conclusion is that the 
level of the predicted background is more than a factor 
of 100 below the observed background. Calibrating on 
the Milky Way, the relative smallness of the obtained re- 
sults for the cosmic gamma-ray background shows that 
the light Dark Matter annihilation hypothesis is by far 
not ruled-out by the current soft gamma-ray extragalac- 
tic background constraint. 

Note that this conclusion can be applied more gener- 
ally to other positron sources since the 511 keV emission 
from the Milky Way is quite weak compared to the ob- 
served background intensity. In order to make this back- 
ground large, one has to assume that other galaxies have 
much higher positron production rates than the Milky 
Way. Indeed, in this case, the Milky Way would not be 
representative from other halos of the same mass. 



102 - 




o 

> 

^ 10-2 - 
10-3 - 

10-* 



10 100 1000 

hi/(keV) 

FIG. 6: Diffuse background spectrum. Crosses: HEAO data; 
stars: COMPTEL observations. The continuous line is the 
sum of the contributions from Seyferts, SNIa and blazars. 
The dot-dot-dot-dashed line and the dashed line represent the 
revised positron contribution to the cosmological gamma-ray 
background, for a S'-wave cross-section with a NFW profile, 
and a P-wave cross-section associated to a Moore profile, re- 
spectively. 

Using the calibration on the Milky Way, we are now 
able to compute the diffuse background for our two 
best models: S'-wave cross-section, NFW profile with a 
100 MeV Dark Matter mass; and P-wave cross-section, 
Moore profile with a 1 MeV Dark Matter mass. These 
are used here as specific benchmarks for the the purposes 
of our analysis, many other intermediate situations being 
obviously also possible. 



IV. CONSTRAINTS ON DARK MATTER 
CANDIDATES 

The main objective of this paper is the computation 
of the Dark Matter induced background in the 100 keV- 
511 keV energy range taking into account the important 
role of baryons. Using standard concentration parame- 
ters for the Dark Matter halo profiles (as given by N body 
simulations), the net result of this exercice is a decrease 
by a factor of 10 of the level of the background emission 
compared to the precedent computation in |l3l |. 

Calibration of {avrd (r)) and on the INTEGRAL 
signal is however uncertain. Indeed, it depends strongly 
on the Milky- Way Dark Matter profile which is not well 
known. The grey region in Fig. [7| shows the range of al- 
lowed parameters corresponding to this calibration. To 
obtain this domain we have considered a reasonable range 
of Dark Matter profiles from peaked and concentrated 
ones (Moore profile and c = 15) to less peaked and less 
concentrated ones (NFW and c = 5). Then we have 
found for each profile which cross-section and particle 
mass reproduce the level of the galactic emission. As a re- 
sult, each mass and cross-section in the grey region could 
reproduce the level of the bulge emission with a reason- 
able Dark Matter profile. The corresponding background 
is below the observed background (see the previous Sec- 
tion). As a consequence, no light Dark Matter candidate 
is ruled out by the diffuse background constraint and one 
has to invoke another gamma-ray source in order to ex- 
plain the missing 4 keV cm~^ s~^ sr"-'^. 

If we consider now that the Milky- Way Dark Matter 
halo is totally different from other halos of similar mass 
in the Universe, we could relax the previous calibration 
and obtain independent constraints on the mass and the 
cross-section (following in that sense the strategy used in 
|l3j | on the observed background only). In the present pa- 
per, constraints are of course less stringent than in |l3j |. 
since baryon physics has led us to decrease the level of the 
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FIG. 7: Constraints on the dark matter candidate in the 
a — mx plane, with a the velocity-independent part in the 
annihilation cross-section a «rei and mx the mass of the dark 
matter particle, h is fixed by the relic density requirement, so 
that {pv^c\)p — < a + bv^ >f ~ 10"^^ cm'^ s~^ at freeze-out. 
The upper line a ~ 10~^^ cm'^ s~^ corresponds to a purely 
S'-wave cross-section, and the lower part of the diagram to a 
P-wave dominated one (a being negligible). As mx cannot 
be too large, in order not to overproduce gamma and radio 
continuum from the galactic center, we limit ourselves, con- 
servatively, to the interval 0.511 to 100 MeV (the actually- 
allowed mass interval could in fact be significantly smaller 
[33 . 145I I4M li^ depending on how seriously these other con- 
straints are taken). 

The grey region is the one compatible with the galactic con- 
straint, based on the total level of emission, using the (Moore 
or NFW) dark-matter distributions of Table I (its upper part 
corresponds to a Milky- Way emission that would-be S'-wave- 
dominated, with a behaving like 1/mx)- 

The dot-dashed line is associated with a dark-matter induced 
background (for a NFW profile) that would correspond to the 
missing 4 keV cm~^ s~^ sr~^, the top left-hand corner above 
this line being excluded on the basis of the cosmic background 
data (same for the dashed line with Moore profiles). 



plained gamma-ray background emission at most of the 
order of 4 keV cm~^ s~^ sr~^ remains. As proposed 
in 11], the strong 511 keV emission from the galac- 
tic center detected by the INTEGRAL satellite could 
be explained by light Dark Matter annihilation, and we 
have verified that the observed emission profile can be 
reproduced, both for S and P-wecve annihilation cross- 
sections. Using the hierarchical model of structure for- 
mation, we have computed the corresponding gamma- 
ray background, and found it to be compatible with 
current observational bounds, ij one takes into account 
the minimal halo mass for galaxy formation. The new 
positron-generated (Dark Matter-induced) extragalactic 
background is in fact overwhelmed by other emissions 
from SNIa, Seyferts, and blazars. The exclusion of small 
mass halos as (redshifted) 511 keV photon sources leads 
to an order of magnitude decrease of the extragalactic 
flux around 500 keV as compared to earlier studies. The 
spectral shape of the extragalactic background is also 
modified in the sense that the number of Dark Matter 
halos capable of hosting gas rich galaxies decrease very 
strongly with increasing redshift. 

Acknowledgements: The authors would like to thank 
the anonymous referee for his helpful remarks that have 
greatly improved the quality of the paper. 



emission by a factor of 10. As shown in Fig. [3 constraints 
exclude only the upper left hand corner, corresponding to 
low-mass Dark Matter particles and S'-wave dominated 
cross-section. Interestingly, for halos with a Moore pro- 
file, the Dark Matter candidate mass and cross-section 
required to reach the level of the observed background 
are not so far from the region favored by the galactic 
constraints. 



V. CONCLUSIONS 



Having estimated the SNIa contribution in the 100 
keV-10 MeV energy range, we have found that an unex- 
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APPENDIX A: 511 KEV EMISSION FROM THE 
GALACTIC BULGE 

The computation of the Milky- Way emission profile 
F{9) is not the main goal of this paper, since we only used 
the integrated emission (and not the shape) for calibra- 
tion purposes. However, it is essential to compare care- 
fully predicted profiles to the observed one (see the recent 
paper '34]). In this Appendix, we would like to outline 
some interesting issues concerning this point. Using the 
same hypothesis as in Section 2, we consider the three 
different cases of Tabled a NFW profile with c = 10 and 
a 5- wave cross-section; a NFW profile with c = 15 and 
a P-wave cross-section; and finally, a Moore profile with 
c = 10 and a P-wave cross-section. 

In the pure S'-wave case with an essentially constant 
((TWrci) —aw 10~^^ cm'^ s~^, we have chosen to con- 
sider and test, with the above NFW distribution, a Dark 
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FIG. 8: INTEGRAL/SPI overall counting rate in the 511 
keV line as a function of the angular distance to the galactic 
center. Points with error bars show the instrumental back- 
ground substracted count rate of the SPI camera in the 506 
to 516 keV band, averaged over the first year of data. 

Lines indicate the expected counting rates for three difi'er- 
ent situations: S-wave annihilation cross-section with a NFW 
Dark Matter profile (continuous line) ; or P-wave cross-section 
with a NFW (dashed line) or Moore (dotted line) profile, as 
specified in the text. The figure does not show the actual 511 
keV emission profiles, but the results of their convolution with 
the spectrometer's response, the Point Spread Function (PSF) 
- which is a crude but more rigourous way in order to com- 
pare data with models. The lower dot-dashed line shows the 
rate (not compatible with the data) that would correspond to 
a pointlike source. For the calibration, we choose the Dark 
Matter particle mass so that the observed integral of the sig- 
nal over the inner 16 degrees is equal to the one from the 
theoretical profile. 



Matter particle with mass mx = 110 MeV/c^. Let us 
note however that, given our hypothesis, the same emis- 
sion profile would have been obtained from the same px 
but with a mass of ~ 11 (or 1.1) MeV/c^ only, and 
an a term that would be ~ 100 (or lO'*) times smaller. 
The corresponding (S'+P-wave) cross-section would then 
be P-wave dominated at freeze-out, while appearing as 
S- wave- dominated for low-velocity annihilation in the 
galactic center. 

In the pure P-wave case on the other hand, for 
which the annihilation cross-section in the galactic cen- 
ter {avrei} should be lower (i.e. typically ~ 10"'^° cm"^ 
s~^), mx should in general be taken relatively small, to 
get (with a correct relic density) a sufficiently intense 
gamma-ray line. In practice we test mx about 0.5 (with 
the NFW profile) and 1.2 MeV/c^ (with the Moore pro- 
file). 

At first sight, all three tested profiles seem compat- 



ible with the observations, within the precision of the 
analysis. Attempting to discriminate between them (or 
with analogous ones) would require a careful chi-squared 
analysis and is beyond the scope of this article, as we 
are mainly interested here in the total observed intensity 
and global morphology of the 511 keV emission of the 
galactic bulge. 

If the cross-section for Dark Matter annihilation in ha- 
los is velocity-independent (S'-wave or effectively S'-wave 
annihilation) , the emissivity of a NFW Dark Matter halo 
scales near the center as P^, (x r~^. Convolving this emis- 
sion with the SPI Point Spread Function (PSF), we ob- 
tain the profile presented as a continuous line in Fig. |S1 
We emphasize here that the profiles shown are the pro- 
files after convolution by the PSF. 

If, however, this cross-section is S'-wave suppressed, 
the emissivity now depends on the Dark Matter 3D ve- 
locity dispersion [s^ . We then compute a^^j as a function 
of the radius by solving the Jeans equation, at first for 
a NFW potential. The resulting emission profile (shown 
after convolution as a dashed line in Fig.|SJ) also turns out 
to fit the data, although it is less peaked, as the velocities 
increase with r, within the region of interest. Further- 
more, as mentioned earlier, to get in this case the appro- 
priate intensity for the 511 keV line we need to consider 
both rather small values of mx (about 0.5 MeV/c^) and 
somewhat extreme parameters for the Milky -Way Dark 
Matter halo (choosing c ~ 15). 

These restrictions may be avoided to some extent, how- 
ever, with a steeper profile such as the Moore profile |30j| . 
which enhances the rate of Dark Matter annihilation, es- 
pecially near the center of the galaxy, so that the result- 
ing emission profile (shown after convolution as a dotted 
line in Fig. IS)) gets now more peaked. With such profiles, 
P-wave annihilation with standard Milky- Way param- 
eters appear to be compatible with the data, even for 
less small values of mx- Furthermore, note that if Dark 
Matter is subject to the stellar gravitational potential that 
dominates the central region of the galaxy (43l , w ith a ra- 
dial density profile declining roug hly as the 5- 
and P-wave cases would be essentially indistinguishable. 

To conclude this Appendix, we have verified that one 
can reproduce the photon flux and distribution observed 
by INTEGRAL, both for S- and P wave cross-sections, 
with standard Dark Matter profiles px and appropriate 
mass mx . Again, in this approach, attempting to further 
discriminate between emission profiles associated with 
S- or P-wave annihilation appears as difficult, given the 
width of the PSF function and the variety of the Dark 
Matter profiles, gravitational potential profiles and gas 
density profiles which may be considered. 



